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Abstract 

There exists a well known problem with the 7 Li+ 7 Be abundance predicted by 
standard big bang nucleosynthesis being larger than the value observed in popula- 
tion II stars. The catalysis of big bang nucleosynthesis by metastable, tx ^ 10 3 sec, 
charged particles X~ is capable of suppressing the primordial 7 Li+ 7 Be abundance 
and making it consistent with the observations. We show that to produce the cor- 
rect abundance, this mechanism of suppression places a requirement on the initial 
abundance of X~ at temperatures of 4 x 10 8 K to be on the order of or larger than 
0.02 per baryon, which is within the natural range of abundances in models with 
metastable electroweak-scale particles. The suppression of 7 Li+ 7 Be is triggered by 
the formation of ( BeX~) compound nuclei, with fast depletion of their abundances 
by catalyzed proton reactions, and in some models by direct capture of X~ on 7 Be. 
The combination of 7 Li+ 7 Be and 6 Li constraints favours the window of lifetimes, 
1000 s £ t x < 2000 s. 



1 Introduction 



Big Bang Nucleosynthesis (BBN) is the process of light element formation during the early 
stages of cosmological expansion. Theoretical predictions of elemental abundances are based 
on known physics such as the nuclear reactions of light elements, general relativity and the 
Standard Model (SM) of particles physics. If the observations of hydrogen, helium and 
lithium abundances in the present are capable of determining their primordial values, the 
consistency of the entire theoretical framework can be tested. In recent years these tests 
acquired particular sharpness, as the only free parameter, the ratio of baryons to photons, 
is now well- measured through the anisotropy of the cosmic microwave background pQ. 

Perhaps the most exciting prospect of studying the primordial abundances is the pos- 
sibility of testing the combination of Standard Model and general relativity, which we call 
SBBN, or standard BBN. To this end it is important to understand how the non-standard 
physics can affect the outcome of nuclear reactions (see e.g. [2 J for a review). Schematically, 
the BBN equations can be represented as 

dY- dY- 

= -H{T)T-^ = J2( T v Y i + T ikiYkYi + ...); Energy of reactants ~ T < MeV 

H{T) = const x N l J 2 — -, where N, = N boson + -N fermion . (1.1) 

In this formula, are the abundances of light elements, T^- are the generalized (positive 
or negative) rates for creation or destruction of element i with participation of j, k..., H(T) 
is the Hubble expansion rate, Mpi is the Planck constant, and is the number of effective 
degrees of freedom comprised of fermionic and bosonic fields. The radiation-domination 
expression for Hit) = l/(2t) is used. There are several ways in which non-standard cos- 
mological and/or particle physics can affect the freeze-out abundances Yi(T — ► 0) of light 
elements. To date the following generic possibilities have been identified: 

1. Timing of the reactions can be changed, for example by models that have additional 
contributions to H(T) (See Ref. [2] and references therein). Such contributions may 
come from the additional thermally excited relativistic degrees of freedom (historically 
often referred to as "extra neutrino" species), or from any other forms of energy that 
contributes to the total energy density at a significant level, such as "tracker" scalar 
fields [3]. 

2. Non-thermal components to nuclear reactions can be introduced by injection of energy 
during or after the BBN pQ. For example, unstable or annihilating heavy particles can 
cause such injection, leading for example to the break-up of He into D or 3 He , or to 
the synthesis of lithium via out-of-equilibrium processes involving energetic secondary 
projectiles. The most important parameters in these types of models are the amount 
of energy injected at t ~ lifetime, and the fraction of energy injected into hadronic 
degrees of freedom. 



1 



3. New catalyzed thermal channels Ty_, can be opened up if the heavy particle physics 
remnants have the ability to bind to nuclei during the BBN. The consequences of 
the catalyzed BBN (CBBN) scenario with charged relics which are present at the 
time of the BBN, and the corresponding catalytic enhancement of 6 Li by several 
orders of magnitude, was explicitly demonstrated in the recent paper [?]. The main 
parameters that determine the strength of catalysis are the initial number densities 
of such particles and their lifetimes. 

Other notable modifications to the standard scenario include an impact from possible 
inhomogeneity of the Universe at the time of BBN [5], and even more contrived options 
such as time dependent couplings [6J. The latter again affects the BBN predictions mostly 
via changing the timing of such processes as the n/p freeze-out and deuterium formation. 

In this paper we assess the ability of CBBN to reduce the amount of primordial 7 Li. This 
question has gained importance recently, as the amount of 7 Li predicted by SBBN is about 
a factor of 2 — 3 larger than the amounts of 7 Li observed in the stellar atmospheres at low 
metallicities. This discrepancy is difficult to attribute to the nuclear rate uncertainties in 
the context of the standard scenario [HI El- It is entirely possible, however, that unaccounted 
depletion of 7 Li in stars can be responsible for the overall depletion of primordial 7 Li. We 
refer the reader to the recent discussion in astrophysical literature [10j . Regardless of 
whether the primordial lithium overabundance issue can be resolved by careful analysis 
of realistic stellar models, one should investigate the feasibility of other options for the 
reduction of primordial lithium, such as by new models of particle physics. It has been 
shown recently that the reduction of 7 Li from a large energy release is possible if the 
metastable particles decay at around 30 keV temperature (tx ~ (1 — 3) x 1000 sec) and 
have a significant hadronic fraction in their decay products [llj . Another possibility based 
purely on the catalysis of thermal nuclear reactions was pointed out in [7|, where a similar 
range of lifetimes has been suggested. Besides Ref. [7j, the cosmological consequences from 
the bound states of nuclei with charged particles were considered in [121 EH1 US] (see 
also the discussion in the earlier work [16J). 

The efficiency of the catalytic reduction of 7 Li, or more precisely of 7 Li+ 7 Be, depends 
crucially on the rate of the ( 7 BeX~) bound state formation , as more than 90% of predicted 
lithium at r/b = r\{f MAP comes in the form of 7 Be. Once ( 7 BeX~) is formed, a new set of 
reaction channels opens up [7], and this paper considers them in detail. These channels 
include an enhancement in p-destruction of 7 Be, enhanced internal capture of X~ and 7 Be, 
as well as other channels. The ( 7 BeX~) bound state then serves as a "bottleneck" for a 
potential solution to lithium problem in the CBBN framework. It is therefore very impor- 
tant to calculate the efficiency of the ( 7 BeX~) formation as accurately as possible, as the 
naive hydrogen-like formula employed in some studies may give a poor approximation to 
the correct answer. Having determined the catalyzed nuclear rates, and the rate of the 
bound state formation, we solve the resulting BBN equations in the lithium chain numer- 
ically for some representative values of the parameter space and determine in which cases 
CBBN provides a noticeable reduction of 7 Be+ 7 Li, while keeping 6 Li within observational 
bounds. 
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It is important to stress that we try to keep our investigation of CBBN processes as 
model independent as possible. The two main parameters that enter in the CBBN reactions 
are the lifetime and the initial abundance of the X~ particles, which we choose to normalize 
on the baryon number density, 

— = Y x exp(-t/r x ) = F x exp(-178s/(r x T 9 2 )). (1.2) 
n b 

In this formula, Tg is the temperature in units of 10 9 K, tx is the lifetime, and Y x is the 
initial abundance at t <C tx- Besides the existence of {Yx, Tx} parameter space, one has to 
recognize the existence of two generic classes of models that can affect the CBBN reactions. 
These are models of type I and type II [7J, which have different ways of achieving the 
longevity of X~. In type I models the long lifetime is due to the weakness of the couplings 
(e.g. models with the gravitino as the lightest supersymmetric particle), while in type II it 
is due to the small, O(MeV), energy release in the transition between X~ and X° while the 
couplings are assumed to be of normal size. Type I models are necessarily accompanied by 
large energy release, and thus possible impact on BBN can be twofold: due to the catalysis 
of thermal nuclear reactions and due to non-thermal processes. In this paper we do not 
treat non-thermal effects, as they tend to be very model-specific. 



2 Recombination of Be and X 

2.1 Properties of the bound states 

The properties of the bound states of 7 Be and X~ are very important because they are 
essential for determining the temperature and probability of recombination. Using the 
Gaussian charge distribution within the 7 Be nucleus, we determine the following energy of 
the ground state relative to 7 Be+X~ continuum: 

Ef ] = -2785 keV; a B = 1.03fm for (r 2 ) Bc = (2.3) 
E g = -1330 ± 20 keV for ((r 2 ) Bc ) 1/2 = 2.50 ± 0.04fm 

In this formula, the first line is a naive Bohr-like formula for a point-like nucleus, which, 
of course, gets a significant correction from the realistic values of the charge radius. These 
values can be derived from experimentally measured values for the nucleon radius of 7 Be 
[17]. In the remainder of the paper we adopt the central value from the second line of (12.31) . 
We have checked that both the variational calculation and the numerical solution to the 
Schrodinger equation produce identical answers. The radial profiles for the wave function 
are plotted in Figure 1. 

Information about the ground state alone is insufficient for calculating the realistic 
recombination rate. The internal structure of the 7 Be nucleus must also be taken into 
account. In particular, the first excited state of beryllium is close to the ground state and 
cannot be neglected, 

|^Be 3/2 - Er B | = 429 keV (2.4) 
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: Radial distance in units of Bohr radius 
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Figure 1: Normalized radial wavefunction for the bound state ( 7 BeX~) as a function of distance in units 
of as = 1-03 fm. A. A hydrogen-like profile for a "point-like" 7 Be. B. Realistic wavefunction for the 
Gaussian charge distribution inside 7 Bc. 

Indeed an excited form of 7 Be, which has total angular momentum I = 1/2 (as opposed to 
the ground state with 1 = 3/2) can also bind to X~ . Some of these metastable bound states 
have total energies just above the 7 Be+X~continuum threshold. This opens the possibility 
of a resonant transition into a bound state, and therefore such states have to be taken into 
account separately. The three most important bound states identified in this way belong 
to the principal quantum number n = 3 

( 7 Be 1/2 X~), n = 3, 1 = 0, E R = (-239 + 429)keV = 190keV; 

n = 3, 1 = 1, E R = (-290 + 429)keV = 140keV; (2.5) 
n = 3, 1 = 2, E R = (-308 + 429)keV = 121keV. 

In these formulae, E R is the energy level of the state relative to the 7 Be 3 / 2 +X~ continuum. 
It is also of interest to quote the position of n = 2 resonances: 

( 7 Be 1/2 X-), n = 2, 1 = 0, E R = (-478 + 429)keV = -49keV; 

n = 2, 1 = 1, E R = (-636 + 429)keV = -206keV. (2.6) 

Unlike the hydrogen-like case, the 21 are not exactly degenerate, which allows for a signifi- 
cant one-photon El transition between them, T2 S ->2p — 0.5eV , and is quite comparable to 
the electromagnetic widths of n = 3 states. One can also notice that the 2s state is very 
close to the threshold, which is potentially very important for the calculation of the capture 
rate. Moreover, there are number of additional effects that may lead to a few tens of keV 
upward shift of this level. Such effects include the correction for the finite mass of X~ , a 
perhaps larger charge radius of the excited state, and correction for nuclear polarizability. 
The precision in our calculations is certainly limited, and only the dedicated many-body 
nuclear studies can reach better than O(50keV) determination of the energy levels. We 
therefore conclude that the 2s state may be indeed at the threshold or just above E — 0. 
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2.2 Resonant and non- resonant recombination 

The process of recombination between 7 Be and X~ may seem to be relatively simple to 
calculate. For the continuum spectrum wave functions we can use the approximation of 
E ~ to good accuracy, as the characteristic energies are on the order of the temperature, 
which is significantly smaller than the binding energy, T/E g ~ 0.03. We find these wave 
functions by solving the corresponding Schrodinger equation numerically. The process of 
recombination can occur via two mechanisms involving resonant and non-resonant capture. 

The non-resonant recombination is given by the processes 

7 Be 3/2 + X- -> ( 7 Be 3/2 X-, n > 1) + 7 -> ( 7 BeX~) + fc 7 , (2.7) 
7 Be 3 / 2 + X- -> ( 7 Be 1/2 X~, n = 1, 2) + 7 -> ( 7 BeX~) + fe 7 , 

where ^ is a typical number of emitted photons ranging from 1 to 3. Among the non- 
resonant processes the capture of 7 Be 3 / 2 directly to the ground state has the largest cross 
section, closely followed by the capture to n = 2, 1 — level. The cross section of pho- 
torecombination differs very significantly from the "Hydrogen-like" formula. In fact, it is 
about one order of magnitude smaller than the naive formula, with most of the suppression 
coming from the u 3 factor , which is a factor of ~8 smaller in the realistic case (12.31) . 

The resonant capture occurs through the following chain of transitions: 

7 Be 3/2 + X- -> ( 7 Be 1/2 X-,n > 3) -> ( 7 Be 1/2 X~, n = 1) + fc 7 -> ( 7 BeX~) + (k + 1) 7 , (2.8) 

To estimate the cross section in this case, one has to find the widths of the corresponding 
metastable states. It turns out that the entrance width of the process 7 Be 3 / 2 + X~ — > 
( 7 Bei/ 2 X~, n = 3) is due to a quadrupole nuclear transition in the electric field of X~. 
If the distance r between X~ particle and 7 Be is large, the interaction takes the form of 
Vi n t = iQi^ViV^a/r). using the existing experimental information on the quadrupole 
matrix elements in 7 Be system, we estimate matrix elements of this interaction between the 
initial wave function of the free 7 Be 3 / 2 , and the intermediate wave function of the 3/ bound 
states of 7 Bei/ 2 , to conclude that Tq is comparable to 1 keV. Given that the electromagnetic 
decay widths of 3/ states, whose calculations are much easier, are all on the order of T° ut ~ 
eV, we arrive at the following hierarchy of the "in", "out" widths and the temperature, 

T out < pin < T ( 2 g) 

This warrants the narrow resonance approximation in the Breit-Wigner formula, and we 
use 

<7R ~ gCTgeom X 7rT^5(E - E R ), (2.10) 

where g is the corresponding multiplicity factor, and a geom is the geometric cross section. 
Notice that in this approximation the rather uncertain value of T 1 q does not enter into the 
formula for the cross section. 
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After performing several quite tedious but straightforward calculations, we arrive at the 
total recombination cross section, averaged over the thermal distribution of X~ and 7 Be: 

6 x 10 3 1.9 x 10 4 . s 1.5 x 10 4 . , , 

(crrecv) = ^- + p— exp(-1.40/r 9 ) + exp (-1 .62/T 9 ) . (2.11) 

^9 ^9 ^9 

Here Tg is temperature in 10 9 Kelvins, and the rate is expressed in conventional astrophys- 
ical units of A^ 1 cm 3 s~ 1 mol _1 . The first term in (12.111) corresponds to the non-resonant 
processes of (12.71) . while the second and third terms are the resonant pieces with internal 
excitations of 7 Be and initial capture on n = 3, / = 2 and n = 3, 1 = 1 levels. The relative 
size of the two effects can be seen in Figure [21 where the total recombination cross section, 
and the non-resonant piece are plotted separately. One can see that the resonant effects 
from 3/ levels may provide a sizable (up to 50%) contribution in the most important domain 
of temperatures, Tg ~ 0.4. 

So far the possible contribution of n = 2 states has been ignored. There is, however, an 
ample chance of 2s state giving a large contribution to the recombination cross section via 
the following chain, 

7 Be 3/2 + X- -> ( 7 Be 1/2 X~, n = 2,1 = 0)^ ( 7 Be 1/2 X-, n = 2, 1 = 1) + 7 -> ( 7 BeX~) + 3 7 . 

(2.12) 

As was argued in the previous subsection, only a dedicated nuclear calculation can exactly 
determine the position of this resonance, which consequently bring a sizable nuclear uncer- 
tainty into the calculation of recombination cross section. In what follows we assume two 
extreme situations. One is when the 2s level is sufficiently below the threshold so that it 
does not contribute significantly to the non-resonant part of the rate, in which case Eq. 
(12.111) will be adopted. Another case, is when we assume that the 2s level is pushed upward 
to the resonant energy of +10 keV. In this case we estimate that the total recombination 
rate from (I2.1ip receives an additional contribution of 

4 x 10 3 

(a rec v) -> (a rec v) H ^— exp(-0.12/T 9 ), (2.13) 

^9 

and the total result of (I2.13P can be viewed as the most "optimistic" estimate of the re- 
combination rate. This contribution is controlled directly by T 2s ->2p- If indeed the resonant 
energy is around 10 keV or so, then for all practical purposes the exponent in (12.131) can 
be taken as ~ 1. As one can see from Fig. (2j the contribution of 2s state is capable of 
enhancing the recombination rate by a factor of a few, thus signaling a significant nuclear 
uncertainty. 

The photodisintegration rate of 7 Be is one-to-one related to the direct rate via the 
detailed balance relation |23j : 

Fphoto = / a photo dn 7 (E) = (a rec v) x I ® e exp(—\E g \/T). (2.14) 

JE>\E g \ V 171 J 
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Figure 2: Recombination rate (cr rec w) for 7 Be and X~ in astrophysical units. A: Nonresonant contri- 
bution. B: Total recombination rate including 3/ resonances. C: total recombination rate, including the 
2s level, in the assumption of Er ~ 10 keV. Given possible O(50) keV uncertainty in the position of 
the 2s level, the whole area between curves A and C is representative of the nuclear uncertainty in the 
recombination rate. 



Numerical expression for the photoionization rate is given by 



^ = {vrecv) x exp(-15.42/T 9 ) , (2.15) 



where n y (T) = 0.24T 3 is the total number density of photons. 

Combining the two rates, (12. lip and (I2.14p . one arrives at the recombination equation 
for the abundance of bound states ( 7 BeX~) relative to the total number of 7 Be nuclei, 

dY 

- HT—^- = n x (l- Y BS )(cr rec v) - T photo Y BS , (2.16) 

where we also assume that n x 3> m Be and took the limit of large t x . One can easily see 
that in the limit T rec ; T photo ^> H this equation has an attractor Saha-type solution, 

Ybs -> Ylf* = [1 + (m Bc T/27r) 3 / 2 n x 1 exp(-|E g |/T)]~ 1 , (2.17) 

which corresponds exactly to the case of chemical equilibrium. 

Solutions for three different values of X~ abundance are plotted in Figure [3] for two 
choices of recombination rate, (12.111) and (12.131) . with and without the contribution of 
the above threshold 2s resonance. It is important to keep in mind that these plots are for 
illustrational purposes only, as they neglect important nuclear effects that destroy ( 7 BeX~), 
which will be considered in subsequent sections. But even this simplified analysis shows 
that the abundance of Y x < 0.005 will be inconsequential for 7 Be abundance as less than 
20% of these nuclei form bound states at T 9 > 0.2. Comparison of two figures exemplifies 
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the nuclear uncertainty, as the efficiency of the recombination rate may be enhanced by 
a factor of a few. Figure [3] also assumes constant abundance Y x , i.e. the decays of X~ 
are neglected. The generalization for a finite lifetime would result in a sharp drop in the 
number of ( 7 BeX~) bound states at times much larger than tx- 




0.2 0.25 0.3 0.35 0.4 0.45 

T 9 




0.2 0.25 0.3 0.35 0.4 0.45 

T 

Figure 3: Fraction of 7 Be locked in the bound state ( 7 BeX~) as a function of temperature. Top figure 
corresponds to a conservative choice of recombination cross section (|2.11[) and to three different values of 
Yx = nx/n-B'- 0.1, 0.03 and 0.01. Lower figure corresponds to recombination cross section enhanced by 
the 2s resonance (|2.13[) and Yx = 0.03, 0.015, 0.005. The long lifetime of X~ is assumed. 



3 Catalyzed reaction channels 

Once the bound state ( 7 BeX~) is formed, several destruction channels get enhanced and 
new destruction channels open up. 
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3.1 Catalysis of p- destruction. 



The reaction Be+p — > 7+ 8 B is one of the best studied in nuclear physics. Its rate is dom- 
inated by an El transition, and consequently depends on the third power of the frequency 
of the emitted photon. As a starting point for our estimates, we use the result of Ref. [18J 
for the value the astrophysical S'-factor at zero energy based on a two-particle model, 

WO) = y« ( ZlA2 - Z2Al ) 2 u? SBBN l\0), (3.18) 

where 1(0) is the radial nuclear integral, and Ai are charges and masses of nuclei 
participating in 7-fusion. 8 B produced in this reaction rapidly undergoes /3 + -decay to 
8 Be thus resulting only in two extra a particles. The SBBN rate for this reaction is not 
sufficiently fast to reduce the amount of 7 Be, as it is about two orders of magnitude slower 
than the Hubble rate. 

The rate for this reaction would have been tremendously enhanced due to virtual photon 
exchange [7] if ( 7 BeX~) + p — > 8 B + X~ were energetically allowed. In reality, this process 
can only go from n > 4 excited states of ( 7 BeX~) and only a small fraction of ( 7 BeX~) 
would be in these states. This does not mean, however, that the rate for the reaction 
involving the ground state is not enhanced. The process 

( 7 BeX-)+j9-^ ( 8 BX-)+ 7 , (3.19) 

produces a photon which is ~ 4.8 times more energetic than in the non-catalyzed reaction, 
and has a reduced Gamow suppression. Assuming that nuclear radial integrals (stripped 
of their Coulomb suppression) have similar values in both the SBBN and CBBN cases, one 
arrives at the following, admittedly crude, estimate for the enhancement of the catalyzed 
S-factor for the reaction in (13.191) : 

Scbbn(O) ~ S S bbn(0) _L ~ 700 x S SB bn(0) ~ 15keV bn, (3.20) 

V^sbbn/ 0.37^ 

where for the SBBN S'-factor we use Ssbbn(O) = 21 eV bn, and 0.37 accounts for the change 
in (ZiA 2 - Z 2 A 1 )/A factor. 

It turns out that in addition to the non-resonant part of the cross section, there are 
more important contributions from resonances. The most important resonant process is 
given by the following transition, 

( 7 BeX") +p -> ( 8 BX", n = 2,1 = 1) -> ( 8 BX~) +7; E R = 167 keV, (3.21) 

where the resonant energy is given relative to the ( 7 BeX~)+p continuum and the input 
value of 2.64 fm is used for the charge radius of the 8 B nucleus. The El electromagnetic 
width of the ( 8 BX~, n — 2, 1 — 1) resonance can be easily calculated to be equal to 
approximately 10 eV. Using this information, and assuming the narrow resonance, we can 
predict the thermal rate for this process as 

(a p v) ~ 1.6 x 10 8 T 9 " 2/3 exp(-8.86/T 9 1/3 ) + 1.6 x 10 6 T 9 " 3/2 exp(-1.94/T 9 ). (3.22) 
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This is a very large rate for a reaction involving 7 in the final state. The resonant part 
of this reaction dominates, and the p-burning of ( 7 BeX~) remains faster than the Hubble 
rate until Tg = 0.2. 

Once the ( 7 BeX~) has being converted to ( 8 BX _ ), the two important processes may 
occur. One is the reverse rate ( 8 BX~) + 7 — > ( 7 BeX~) + p, that brings back 7 Be, and 
the other process is the /5-decay of ( 8 BX~), for which we will use the standard lifetime 
value of 8 B, tb = 1.11 sec. At high temperatures, the reverse rate reduces the efficiency of 
p-destruction very significantly, as pointed out by K. Jedamzik in larXiv:0707.2 070. Since 
the destruction of ( 8 BX~) is very rapid compared to the Hubble scale, either by the beta 
decay or by the reverse reaction T 7 , one can employ the dinamical equilibrium assumption 
for the concentration of ( 8 BX~), which leads to the following effective destruction rate of 
( 7 BeX~) by the proton reactions, 

Tf (T) = (3.23) 

Calculating T 7 according to the standard procedure [23J, we plot the effective destruction 
rate in units of Hubble rate in Figure HI 

Tf/H 

6: 
5: 
4: 

3 = 
2 = 
1 = 
= 

0.1 0.15 0.2 0.25 0.3 0.35 0.4 

T 9 

Figure 4: Effective rate for proton destruction of ( 7 BcX ), that becomes efficient if the ratio T° s /H is 
on the order or larger than unity. Even though the rate for ( 7 BcX~) + p — > ( 8 BX~) + 7 is much faster 
than Hubble rate, the reverse rate makes this destruction mechanism inefficient for T > 30kcV. 

Figure H] shows that in interval of temperatures from 20 to 30 keV the proton burning 
of ( 7 BeX~) is very efficient. After (3-decaj most of ( 8 BX~) bound states would transform 
to ( 8 BeX - ), which would remain stable until the decay of X~, which in turn will be 
immediately followed by the decay 8 Be— *-2 4 He. As an interesting side remark, we observe 
that intermediate ( 8 BeX~) bound states could potentially lead to a new primordial source 
of 9 Be via the catalyzed ( 8 BeX~) + n — > 9 Be + X~ reaction. 
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3.2 Catalysis of n-dest ruction 



It is well-known that the reaction 7 Be + n —> 7 Li + p is the main mechanism for depleting 
7 Be abundance in SBBN. This is a reaction with one of the largest cross sections known 
to nuclear physics with A < 8. It is dominated by a wide near-threshold resonance in 
8 Be (E»Be = 18.9MeV) with the neutron and proton widths being 0.22 and 1.41 MeV 
respectively |19j . With such large widths, the binding effects to X~ are not going to 
change the rate for the ( 7 BeX~) + n — > ( 7 LiX _ ) + p reaction in any substantial way. 

The enhancement of the 7 Be + n — » 2 4 He was mentioned in [8] as one of the important 
destruction mechanisms if its rate is scaled up by an arbitrary enhancement factor of a 
hundred. Ref. [8] concludes that such rescaling is unlikely in this channel. One of the 
papers that investigated the effects of the bound states with charged particles [13] argues 
that the rate for the ( 7 BeX~) + n —* X~ + 2 4 He process, 

2 x 10 4 (l + 3.7 10 3 T 9 ), (3.24) 

could indeed be greatly enhanced by up to three orders of magnitude because 7 Be on the 
orbit of the ( 7 BeX~) bound state has large kinetic energy that effectively replaces T$ in 
the formula for the rate. The absence of T g in the rate above would necessarily mean s- 
wave rather than p-wave scattering, which can only happen if En = 18.9 MeV resonance 
dominates the cross section. The reason why the s-wave part is relatively small in the SBBN 
rate (13.241) is because the symmetry of the final state forbids the photonless decay to 2 4 He 
and allows only 7 Be + n — ► 8 Be* — > 7 + 2 4 He. It is clear that in the CBBN process, internal 
conversion may occur, enabling a photonless s-wave reaction ( 7 BeX~) + n — > X~ + 2 4 He. 
The overall enhancement of the s-wave (Tg-independent) part of (13.241) is expected to be 
on the order of (A rca i/A virtua i) 3 , which for this transition is not likely to exceed 10 3 . Based 
on this argument, we conclude that the rate for ( BeX~) + n — > X~ + 2 4 He is catalyzed 
rather moderately, and would not exceed ~ fewxlO 8 . Such a rate is still subdominant to 
7 Be + n — > 7 Li + p, and the destruction of ( 7 BeX~) by neutrons is less efficient than by 
protons (I3.22p . especially taking into account that neutrons are less abundant than protons 
by approximately six orders of magnitude at relevant temperatures. 

3.3 Break-up of 7 Be in the decay of X . 

New channels of 7 Be break-up come from the decay of the X~ particle inside the ( 7 BeX _ ) 

bound states. There are three different possibilities: 

A: ( 7 BeX~)^ 7 Be + products of decay + hard 7 — > 3 He+ 4 He+... 

B: ( 7 BeX-)-^ 7 Be(£ > T) + background particles + ...-> 3 He+ 4 He+... 

C: ( 7 Be*X-)^ 3 He+ 4 He+... 

In case A, the process of decay of X~ and/or the charged products in the decay chain 
produce a photon with energy in excess of the photodisintegration threshold for 7 Be, E 1 > 
E t hr = 1-59 MeV, that hits the 7 Be nucleus and destroys it. This is a hard photon with 
respect to the thermal bath, T ~ O(30 keV) , but is considered a soft photon with respect 



11 



to the energies of the X~ decay in type I models, where the typical energy release is 
0(100 GeV). The process B is much simpler. The decay of X~ occurs instantaneously, 
so that the emergent 7 Be nucleus has a distribution over the range of momenta given 
by the Fourier transform of the wave function from Figure [U Since the typical kinetic 
energies of 7 Be on orbit inside ( 7 BeX~) are on the order of a few MeV, the break-up of the 
recoiling 7 Be by other nuclei in the primordial plasma becomes more efficient as the Gamow 
suppression is effectively lowered. The final option C is due to the polarization of 7 Be by 
the electric field exerted by X~ inside ( 7 BeX~). Effectively, the 7 Be within the bound state 
is represented by the ground state of isolated 7 Be admixed with a combination of all the 
excited states, most of which (with the exception of 7 Bei/ 2 ) belong to the continuum. Upon 
the instantaneous decay of X~ these extra "polarized" pieces of the 7 Be wave function lead 
to a decay to the 3 He+ 4 He continuum. Mechanism C exists for both types of models, with 
or without large energy injection. 

To estimate the efficiency of channel A, we take a model of type I with the assumption 
that the decay produces a single charged track with energy E max ^> 1.59 MeV. We estimate 
the number of Weizsacker- Williams photons dn 1 {E) with E -C E max , determine the effective 
flux for the 7 Be target and obtain the following probability for the radiative break-up, 

POO l-Emax 

P radbr ~ dr#(r)| 2 / a 7 (E)dn y (E) (3.25) 

J0 J E thr 

Numerical evaluation of (I3.25P with the input of cross section for 7 Be + 7 — » 4 He+ 3 He 
produces the following estimate, 

Pradbr-W- 4 for E max = 100 GeV, (3.26) 

which forces us to conclude that the radiative breakup of the recoiling 7 Be is rather ineffi- 
cient. 

Estimates of process B depend very sensitively on the rate of slow down of 7 Be. For the 
relevant temperature range (on the order of 30keV), we determine that it takes approxi- 
mately 10 _6 sec for an MeV-energy 7 Be to thermalize. This time interval is not sufficient for 
a proton induced break-up to happen, 7 Be+p ^ 9 B+7, even though the rate for the latter 
is enhanced relative to thermal rate due to a larger center of mass energy of 7 Be + p system. 
The estimated probability for the break up in mechanism B during the thermalization time 
remains negligibly small, and we conclude that process B is also inefficient. 

Finally, in option C we need to estimate the admixture of excited states to the ground 
state of 7 Be nucleus by the electric field of X~ particle. Schematically, this probability is 

2 

(3.27) 

and the sum runs over all continuum states of 7 Be. Taking a typical value of (dE) i ~ 1 
MeV, and energy denominator comparable to Gamow energy, Ei — En > 10 MeV, we 
estimate that the total probability of admixture of excited states is only about 1%. Thus, 
we are forced to conclude that all mechanisms related to the decay of X~ within the bound 
state ( 7 BeV~) do not lead to a significant depletion of 7 Be abundance. 



(dE) 



0/ 



En — Ei 
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3.4 Internal conversion: Be + X — ► Li + X . 

Any process that converts 7 Be to 7 Li early enough would lead to a suppression of the 
total lithium abundance, 7 Li+ 7 Be, as 7 Li is significantly more fragile. In type II models 
of charged relics, where their longevity is ensured by the small energy splitting between 
X~ and X°, it is natural to expect the existence of weak charged current between X° and 
X~, X + . In that case the process shown in Figure [5] becomes possible, 

( 7 BeX~) -> 7 Li +X°. (3.28) 

This process is analogous to a usual electron capture process, ( 7 Be e~) — > 7 Li + z/. The main 




Figure 5: Internal conversion of beryllium into lithium in type II models due to PF-exchange. 7 Li, being 
intrinsically more fragile than 7 Be , is subsequently destroyed. 

difference is that the rate for (I3.28|) is many orders of magnitude larger than the standard 
capture of e~ inside the beryllium atom. There are two primary reasons for that. The 
usual capture rate is proportional to 1^(0) | 2 ~ (ZB e am e ) 3 , i.e. to the probability density of 
finding an electron at the same point as 7 Be. For the reaction (13.281) this factor should be 
taken to be ~ , since the average distance between 7 Be and X~ is less than the nuclear 
radius r Be . The second source of enhancement is in the final state momentum of particles 
involved in the transition. More specifically, due to different kinematics, the square of the 
neutrino momentum p 2 % must be substituted for muPu- Altogether this amounts to an 
extremely large enhancement factor, 

N enh ~ (r Be Z Be am e ) 3 x (m u pu/p 2 {u) ) ~ 10 15 ; T int conv ~ 10" 6 eV. (3.29) 

This must be compared to the rate of the normal capture, which is in the range of 
CKKTSec- 1 )- 0(l(T 21 eV). 

Comparing rate (13.291) with the photodissociation rate of ( 7 BeX~) we conclude that 
(13.291) becomes dominant for T 9 < 1, which covers the whole range of temperatures rel- 
evant for the lithium abundance problem. Because of this, in type II models the rate 
of recombination determines the 7 Be— > 7 Li conversion rate induced by X~ , and the actual 
width (13.291) does not enter Boltzmann equations, as it is larger than the rest of the reaction 
rates. 
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3.5 Energy injection through the annihilation of X and X + . 

The impact of annihilating particles on the BBN predictions was recently considered in 
detail in Ref. [20]. In the CBBN scenario, the injection of energy due to the annihilation 
of X~ and X + can be calculated with reasonable accuracy. At temperatures much below 
the binding energy of the (X~X + ) system, this annihilation proceeds via the formation 
of a positronium-like bound state. This recombination process is followed by the rapid 
annihilation of X~ and X + within the bound state. The value of binding energy = 
mx<y 2 /4: is larger than 1 MeV for nix > 100 GeV. Therefore, at T < 100 keV it is safe to 
ignore the reverse process of photodisintegration of (X~X + ). 

A straightforward calculation gives the amount of energy released per X~ particle per 
time, 

Tann = {<T{X~X+)V rel)n X + (3.30) 

Unlike the case of 7 Be+X~ recombination, the hydrogen-like (or more precisely, the positronium- 
like) expressions for the cross sections are appropriate in this case, 

2 9 vr 2 a 2El I 1 \ 2 10 7r 3 / 2 a 3 
<"C*-* + >*"> " 3^)^7 W " 3exp(4) m fTV 2 ' (3 " 31) 

where /x is the reduced mass, /i = mj/2, and Ej, = a 2 mx/^- Notice that it is the mass 
of the X-particle, rather than the nuclear mass, that sets the scale of recombination cross 
section, and consequently the rate f 1 3 . 3 j) is subdominant to the Hubble rate. Therefore the 
formation of (X + X~) positronium-like state with its subsequent annihilation does not lead 
to an appreciable depletion of nx- It is instructive to calculate the actual size of this rate 
at temperatures relevant for the lithium problem: 

(500 Gtg"V \ ^/^ 
— rn~ — J ' ( ' 3 ' 32 ' ) 

This is to be compared with the sub-picobarn size of cross section for the annihilation of 
neutralino particles. Therefore the effect from (X + X~)-induced energy injection is expected 
to be much more important for the BBN than neutralino annihilation [20] . 

It is customary to quantify the effects of the unstable particles on BBN by the energy 
released per photon £ normalized on 1 GeV. A similar quantity for annihilating particles 
would read as 



2m X v2 2 f 2 WiX'X+^rel) 



where T\ to T 2 is the relevant temperature interval, which for our estimates can be taken 
to be from ~ 10 to 40 keV. Above these temperatures, the impact of energy injection will 
be washed away by thermal nuclear reactions. 

Substituting (I3.3ip into (13. 33ft . and normalizing Yx and mx on their typical "electroweak 
scale-inspired" values, we arrive at 

12 /500 GeV\ 1/2 / Y x \ 2 
£ = 2.2 x 10~ 12 — - . 3.34 
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This small energy release per photon is inconsequential for the BBN predictions if the energy 
is released via the electromagnetic radiation, e.g. X~+X + — >(X~X + ) — > 77. However the 
hadronic energy release via X~+X + — *(X~X + ) — > qq even at this O(10 -12 ) level is capable 
of suppressing 7 Li. For example, provided that the hadronic branching ratio is order 1, the 
estimate H3.34I) is consistent with a factor of a few decrease in 7 Li abundance reported in 
[TTj . It is important to note that the realistic scaling of (13.341) with mass is m 3 ^ 2 , as the 
standards arguments for the annihilation at the freeze-out force Yx to scale linearly with 
m>x when it becomes heavier than the electroweak scale. Thus for heavier mx the effect is 
even more pronounced. 

Subsequent decays of X~and X + may or may not lead to a larger energy release than 
residual annihilation. In models of type I, the decay would typically provide larger effect, 
while in models of type II the annihilation is more important. In any event, this subsection 
shows that even for type II models, in which the decays of X~ and X + do not release a 
significant amount of energy, the non-thermal component of BBN is still unavoidable if the 
hadronic branching in the annihilation of (X~X + ) is appreciably large. 

4 Evolution of 7 Be + 7 Li at T ~ 30 KeV 

The standard mechanism of generating 7 Be+ 7 Li at the WMAP-suggested value for r)b is 
relatively simple. The main reactions that determine total lithium abundance are 

SBBN : 3 He + 4 He -> 7 Be + 7; 7 Be + n -> 7 Li + p; 7 Li + p^ 2 4 He or D + 6 Li. (4.35) 

The first reaction generates Be while the combination of the second and third burn it. The 
secondary creation mechanism is due to 4 He+ 3 H— > 7 Li+7, while the secondary destruction 
mechanisms are given by 7 Be + p — > 8 B + 7 and 7 Be + D — > p + 2 4 He. All of these processes 
were discussed at length in the existing BBN literature [H [22] . 

The CBBN creates several additional destruction mechanisms that were described in 
the previous sections. Below we list main CBBN processes: 

CBBN: 7 Be + X" <-> ( 7 BeX~) +7; 7 Li + X~ <-> ( 7 LiX") + 7; (4.36) 

Type I and II : ( 7 BeX") + p <-> ( 8 BX~) + 7; ( 8 BX") -> ( 8 BeX"). 

( 7 BeX~) + n -> ( 7 LiX~) + p; ( 7 LiX") + p -> X' + 2 4 He or X~ + D + 6 Li. 

Type II only: ( 7 BeX") -> 7 Li + X°; ( 8 BX") -> 8 Be + X°. 

Besides these main channels, we also included the ( 7 BeX~) destruction by D-burning. In 
( 7 LiX~) +p and ( 7 BeX _ )+D reactions the only change implemented relative to the SBBN 
rate was in the Coulomb penetration factor. 

The models of Type II have an additional advantage over type I because of the internal 
conversion process ( 7 BeX~) — > 7 Li + X°. As soon as the rate for this process becomes 
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dominant over the Hubble rate, one can consider the photorecombination of 7 Be and X 
as leading directly to 7 Li, 

7 Be + X- -> ( 7 BeX-) + 7 -> 7 Li + X° + 7, (4.37) 

with the rate given by Eqs. (12.111) or (I2.13p . 

The actual calculations of the abundance are performed using only the Li-Be part of 
the BBN code. Since the abundances of these elements are extremely tiny, their "backre- 
action" on 4 He, D and 3 He is negligible. That is, the destruction or creation mechanisms 
of 7 Be and 7 Li do not manifest themselves in any noticeable change in the lighter elemental 
abundances. Thus the CBBN effect on lithium and beryllium can be studied separately, 
and we do so by creating a corresponding code with the use of MAPLE 9.5. The network 
of reactions given in (14.351) . which we take from well-established sources |23j, and the input 
of D(T), 3 He(T), etc. from the full SBBN code gives an approximation to the freeze-out 
abundance of total lithium as 

( 7 Li tot ) S BBN = ( 7 Be + 7 Li) S BBN = 4.1 x 1(T 10 , (4.38) 

in excess of the Spite plateau value, and indeed in agreement with more elaborate treatments 
ofRefs. 0EUE2]. 

We then extend our code by inclusion of the CBBN reaction chains, and do so separately 
for Type I and Type II models. In this paper we do not treat the non-thermal BBN 
processes due to the energy release associated with the decay and annihilation of X~ and X + 
particles. Such processes indeed could be important in type I models, and lead to significant 
modifications of elemental abundances. They are inevitably more model dependent, and 
were previously considered elsewhere (HEED [15]. We also perform calculations for two forms 
of the recombination rate, with and without the effect from the intermediate 2s state in 
( 7 BeX-). 

The temperature evolution of 7 Be, 7 Li and their bound states with X~ are shown for 
typical values of the lifetime and abundance of X~ is presented in Figures [6] and UJ The 
total freeze-out abundance of lithium is changed in the following manner: 

Eq. (2.11); Y x = 0.05; t x = 2000s 

Eq. (2.13); Y x = 0.05; r x = 2000s 

where (14.391) corresponds to a conservative choice of recombination cross section ( 12.111) . 
while (I4.40p takes into account the contribution of 2s resonance H2. 13f) in the assumption 
that this resonance is just above the energy threshold. The overall reduction of lithium 
abundance for this choice of parameters is 10% in type I model and almost 50% for the 
type II for the conservative value of recombination cross section, and 25% and a factor of 
~ 3 for the recombination cross section enhanced by the 2s resonance. 



7T .tot _ / 3.7 xlO- 10 for type I 
LlcBBN _ \ 2.3 x 10- 10 for type II l4 '^ j 

7 tot _ / 3.3 x 10~ 10 for type I , , 
LlcBBN _ 1 1.5 x 10- 10 for type II l4 ' 4Uj 
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Figure 6: The abundance of lithium and beryllium in two types of the CBBN models for the sample 
values of the CBBN input parameters, Yx = 0.05 and tx — 2000 sec, and the choice of recombination cross 
section (|2~TTj) . The total lithium abundance is 3.7 x 10~ 10 for Type I and 2.3 x 10~ 10 for Type II model. 
The individual abundances of 7 Be, 7 Li, ( 7 BcX~) and the SBBN curve for 7 Li+ 7 Be are also shown. All 
abundances are given relative to hydrogen. 
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.1 .2 T 9 -4 

Figure 7: The abundance of lithium and beryllium in two types of the CBBN models for the sample 
values of the CBBN input parameters, Yx = 0.05 and tx = 2000 sec, and the choice of recombination cross 
section (|2"Tl3]l . The total lithium abundance is 3.3 x 10~ 10 for Type I and 1.5 x 10~ 10 for Type II model. 
The individual abundances of 7 Be, 7 Li, ( 7 BeX~) and the SBBN curve for 7 Li+ 7 Be are also shown. All 
abundances are given relative to hydrogen. 
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The behavior of individual curves in Figures [6] and [7J can be easily explained. The 
SBBN curve remains almost flat below Tg = 0.4, while the CBBN introduces an extra 
suppression, resulting in a slight decrease of the total lithium abundance below Tg = 0.4. It 
is important to note that the overall abundance of ( 7 BeX~) remains very low in models of 
type II, being controled by the fast internal capture process. In type I it can be significant 
at 0.25 < Tg < 0.4 but then the p-burning of ( 7 BeX~) makes a noticeable "dent" in 
the ( 7 BeX~) abundance around T 9 ~ 0.25. At later time, all X-containing nuclei decline 
exponentially due to the decay of ( 7 BeX~) at T 9 < 0.2. At the same time, 7 Li abundance 
increases slightly below Tg = 0.3 because the p-burning of 7 Li becomes less efficient. 

Next we explore the parameter space of possible {Y x ,t x } values that would suppress 
total lithium abundance to an "acceptable" level. It is clearly a rather contestable issue of 
what to call an "acceptable level". We take this value to be 2.5 x 10~ 10 leaving some room 
for further possible reduction of lithium abundance in stars. The results of the scan through 
the parameter space is given in Figure [HJ This figure presents to models, with and without 
internal conversion, separately. One can see that the Yx ~ O(.05) and Tx ~ (1 — 10) x 10 3 s 
can provide an appreciable reduction of total lithium abundance. However the catalytic 
production of 6 Li at Tg ~ 0.09 reduces the available parameter space quite significantly. 
Approximating the 6 Li bound [7J by 

Y X (T 9 = 0.09) < 10~ 6 =>■ Y x x exp(-2.2 x 10 4 s/r x ) < KT 6 , (4.41) 
we find that only rather short-lived X~ particles, 

few x 100 s < t x < 2000 s (4.42) 

can provide a reduction of 7 Be+ 7 Li abundance without producing too much 6 Li via the 
formation of ( 4 HeX~) bound states. Overall, once can see that the inclusion of internal 
conversion of ( 7 BeX~) to 7 Li and X° increases the efficiency of CBBN by a lot, allowing for 
a much more parameter space where both lithium problems are resolved by the catalytic 
means. The figures also show that the inclusion of the 2s resonance can increase the 
efficiency of the CBBN by a lot, and therefore a dedicated calculation of the nuclear levels 
in ( 7 BeX~) system with better than ~ 50 keV accuracy is certainly warranted. For models 
of type I one can see the dramatic decrease in the efficiency of CBBN reduction of 7 Li 
abundance for small values of Tx, while in models of type II this effect is less pronounced. 
It can be directly attributed to the fact that p-buring of ( 7 BeX~) is only efficient for 
Tg ~ 0.25, and most of the X~ with lifetimes of the few hundred seconds would decay by 
then. 

5 Discussion and conclusions 

Thermal catalysis of nuclear reactions by metastable heavy charged particles is a new way 
in which particle physics can affect primordial abundances. In the previous paper [Tj a new 
strong source of 6 Li was pointed out, originating from the recombination of 4 He nuclei and 



19 



-10 




Figure 8: {Yx, tx} parameter space with CBBN Lithium abundance. The top panel corresponds to type 
I models, the lower panel corresponds to type II. Curves A and B are calculated with the use of (12.111) 
and (|2.13p correspondingly. The Be+ Li overabundance problem is solved above curves A and B, with 
( 7 Li tot )cBBN < 2.5 x 10~ 10 . The domain to the right from the almost vertical line is excluded from the 6 Li 
overproduction at Tg = 0.09. 
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X~ at temperatures of about 8 keV. The bound on 6 Li is so strong that it requires all 
X~ particles decay before that, r x < 5000 s. In this paper we analyzed earlier stage of 
nucleosynthesis at about 30 keV temperatures. During that period, the only stable bound 
state that can be formed is ( 7 BeX~). Its formation reduces the total 7 Be+ 7 Li abundance 
because this state is more fragile than 7 Be. It is therefore tempting to speculate that the 
CBBN may provide a satisfactory solution to the lithium overproduction problem using the 
catalysis of thermal nuclear reactions. 

We find that indeed a CBBN reduction of the total lithium abundance is possible for 
some values of the abundance-lifetime {Y x ,r x } parameter space. This reduction happens 
more efficiently in models with nearly degenerate X~ ,X + and X° states, connected by the 
weak current, enabling fast internal capture process of Fig. [5j For models without the 
internal capture, the main mechanism for 7 Be depletion is the catalysis of p-burning of 
( 7 BeX~) bound states that at 0.2 < T 9 < 0.3 is faster than the Hubble rate. At the same 
time, we note that the problem of calculating 7 Be+ 7 Li abundance in CBBN is in certain 
ways harder than the previously considered case of 6 Li, simply because it requires better 
than 30% accuracy for the answer, while a factor of a few error in 6 Li can still be tolerated 
given the uncertain observational status of 6 Li. 

The allowed part of the parameter space requires r x to be less than about 2000 seconds, 
and the initial abundance to be of order 0.05 or larger. Interestingly enough, a very similar 
range of lifetimes was suggeted in [TTJ where the lithium overproduction problem is solved 
via the hadronic energy release. Even though such abundances can be achieved in particle 
physics models, it is not difficult to see that they require the total energy density carried 
by X~ be equal or larger than the energy density of cold dark matter. Assuming that the 
minimal allowed mass for X~ is on the order of 100 GeV, we get the following estimate for 
the energy density of X~ before their decay 

fix- iymy 

-J- ~ x x > 5. (5.43) 
\\ m p 

Since it is natural to expect that the same amount of energy would be concentrated in X + , 
we come to the conclusion that in order to solve lithium problem X~ and X + would have 
to carry a factor of two more energy than the cold dark matter. This forces a conclusion 
that the decays of charged X particles in these models would have to be accompanied by 
a significant energy release. There is, however, a caveat of the nuclear uncertainty in the 
recombination rate of 7 Be and X~ , and only more elborate nuclear calculations can clarify 
this picture. 

It is important to realize that the efficiency of the CBBN reduction of 7 Li+ 7 Be is one- 
to-one related to the photorecombination of ( 7 BeX~), or in other words ( 7 BeX~) serves as 
a "bottleneck" for the CBBN reduction of lithium, much like D formation is the bottleneck 
for the helium formation in SBBN. The calculation of the non-resonant contribution of 



1 Recent three-body calculation [55] revises the original estimate of the catalytic enhancement of the S- 
factor from about 1.6 x 10 7 to 2.1 x 10 6 , which does not alter main conclusions of Ref. 7 in any substantial 
way. 
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7 Be capture directly to the ground state of ( 7 BeX~) does not requires significant nuclear 
physics input, and we believe that the current paper provides a reasonable, perhaps ~ 30%, 
accuracy for that part of the cross section. The impact of internal excitations of 7 Be in the 
process of capture may lead to the additional resonant contributions. In particular, we find 
that the 2s atomic resonance between the spin 1/2 (excited) state of 7 Be and X~ is very 
close to the threshold, and may lead to the factor of a few enhancement in the capture rate. 
Unfortunately, we are unable at this point to improve the calculation to the point when 
the contribution of this resonance is firmly established. In practice, one needs to calculate 
the binding energies in the excited states with ~ lOkeV accuracy, which is possible only 
with the advanced techniques of modern nuclear physics. Therefore, such calculation would 
perhaps be the most important rate calculation for the whole CBBN paradigm. 

There is one alternative CBBN scenario that we would like to briefly mention here. 
Suppose there is a pair of closely degenerate massive X particles and stable dark matter 
X° particles that can form a vertex together with an electron line, Fig. [9j This can happen 
if are fermions and X° are bosons or vice versa. If the mass difference of the X~ and 
X° satisfies the following relation, 

m e < m x ~ — %« < 1.33MeV — m e (5.44) 

the two processes are allowed: the decay of X~ and X + to X°, and X°— 7 Be recombination, 

X ± -> X° + e ± ; X° + 7 Be -> ( 7 BeX~) + e + (5.45) 

In this case, the catalytic destruction of 7 Be will proceed the same way as in type I models. 



X 




Figure 9: If the splitting between X~ and X° is between 0.511 and 0.822 MeV, both the decay of X~ 
to X° and the X° capture by 7 Be become possible. At the same time the ( 4 HeX~) bound state cannot 
be formed from X° and 4 He and the 6 Li overproduction is automatically avoided if the lifetime of X° is 
arbitrarily short. 

Since the coupling between e, X~ and X° is in principle a free parameter, one could tune 
it to sufficiently large values so that the second process in (I5.45P exceeds the rate in (12.111) . 
thus enabling more efficient catalytic suppression of 7 Li+ 7 Be. This model would predict 
dark matter capture by nuclei continuing now, which is an intriguing option to consider. 
Detailed analysis of this model falls outside the scope of this paper. 

To conclude, this paper provides the detailed account of main processes responsible for 
the catalytic suppression of 7 Li+ 7 Be in CBBN. The efficiency of such suppression depends 
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mostly on the recombination rate leading to ( 7 BeX~), that is calculated in detail. This 
paper treats only the change in the thermal processes, although the non-thermal effects due 
to the decay and annihilation of charged relic particles are also expected to be important. 
Thus the question of whether or not the lithium abundance can be sufficiently depleted by 
the combination of CBBN processes and the energy release should be addressed on model- 
to-model basis. The stau-NLSP / gravitino-LSP model has already been extensively studied 
in the literature, [26] and [27J. Subsequent work [21] will address the issue of whether the 
significant depletion of 7 Li+ 7 Be can be achieved in the constrained MSSM models with 
stau-NLSP/neutralino-LSP scenario, when the decay of staus is delayed due to the small 
mass difference with neutralino. 
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